Suppression of the rate of growth of dynamic heterogeneities and its 
relation to the local structure in a supercooled polydisperse liquid 
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The relationship between the microscopic arrangement of molecules in a supercooled liquid and 
its slow dynamics at low temperature near glass transition is studied by Molecular Dynamics 
(MD) simulations. A Lennard- Jones liquid with polydispersity in size and mass of constituent 
particles is chosen as the model system. Our studies reveal that the local structure (that varies 
with polydispersity) plays a crucial role both in the slowing down of dynamics and in the growth 
of the dynamic heterogeneities, besides determining the glass forming ability (GFA) of the system. 
Increasing polydispersity at fixed volume fraction is found to suppress the rate of growth of dynamic 
correlations, as detected by the growth in the peak of the non-linear density response function, 
Xii't)- The growth in dynamical correlation is manifested in a stronger than usual breakdown of 
Stokes-Einstein relation at lower polydispersity at low temperatures and also leads to a decrease in 
the fragility of the system with polydispersity. We show that the suppression of the rate of growth 
of the dynamic heterogeneity can be attributed to the loss of structural correlations (as measured 
by the structure factor and the local bond orientational order) with polydispersity. While a critical 
polydispersity is required to avoid crystallization, we find that further increase in polydispersity 
lowers the glass forming ability. 

PACS: 64.70.pm, 61.20.Lc, 82.70.Dd 



I. INTRODUCTION 

The relation between the local structure and 
its slow dynamics in a supercooled liquid near 
glass transition temperature, Tg is currently a 
subject of intense curiosity. The most distinc- 
tive feature of glass formation is the rapid in- 
crease of viscosity with decrease in tempera- 
ture. The temperature at which the viscosity 
becomes 10^'^ Poise is defined as the glass tran- 
sition temperature. One of the main difficul- 
ties in understanding the glass transition phe- 
nomenon is that this enormous slowing down 
of dynamics is apparently not accompanied by 
a growing static correlation length (unlike the 
usual critical phenomena). Static structural 
quantities do not reveal any long range corre- 
lation. In fact the static structure of the liquid 
near glass transition is not much different from 
its equilibrium high temperature coimterpart. 

In the Adam-Gibbs picture the sharp 
slowing down is related to the growth of a co- 
operative dynamic length scale. In a separate 
theoretical study the size of heterogeneous re- 
configuring regions in a deeply viscous liquid 
was inferred via the Random First Order Tran- 
sition Theory (RF0T)[2]. There is now increas- 
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ing evidence from both experiments and simula- 
tions of a dynamic correlation length that grows 
upon approaching the glass transition Q. 
Multipoint susceptibilities have been devised to 
quantify the behavior and magnitude of grow- 
ing dynamic length scales and have been used 
in the experimental studies for several mate- 
rials These have directly determined the 
number of molecular units that move coopera- 
tively near glass transition. The simplest den- 
sity correlation function that contains informa- 
tion on correlated motion is the fourth-order 
0| ■ The four-point time-dependent density 
correlation function, gi{r,t) measure the spa- 
tial correlations between the local liquid den- 
sity at two points in space, each at two different 
times. The dynamical four-point susceptibility, 
Xi{i) (the volume integral of gi{r,t)) becomes 
increasingly pronounced as glass transition is 
approached. 

In this study, we look for a possible relation- 
ship between the structure and the slowdown 
of dynamics in supercooled polydisperse liquids 
near glass transition. In particular, we look at 
how the local structure (which we characterize 
using structure factor and bond orientational 
order parameters) would influence the growth 
of dynamic heterogeneity and the glass form- 
ing ability of the system. Polydisperse liquids 
are one of the simplest model systems that ex- 
hibit glass transition and can be conveniently 
studied via both experiments 0, [1] and com- 



puter simulations as the size distribution of 
particles prevents crystallization [l^, It 
also serves as a model for colloids and many 
other real world systems like polymers, pig- 
ments, paints etc as polydispersity is inherent 
in all these systems. Polydispersity introduces 
a distribution of particle diameters and masses 
and thus makes the system intrinsically more 
heterogeneous. However, the effect of polydis- 
persity on dynamic heterogeneity has not yet 
been examined in detail. Here we probe this is 
detail using the dynamical four-point suscepti- 
bility, Xiit)- Increasing polydispersity results 
in the loss of structural order. Thus by vary- 
ing polydispersity one can understand the effect 
of loss of structure on the growth of dynamic 
heterogeneities. Our studies [ll[ have shown 
that increasing polydispersity at fixed volume 
fraction decreases the fragility. And hence this 
study also presents us with an opportunity to 
probe the growth of four-point susceptibility 
(and thus dynamic heterogeneity) in systems 
with varying degree of fragility. 

The rest of the paper is organized as follows. 
In section // we describe the model and simu- 
lation details and also define various quantities 
that are used in the analysis. In section /// we 
present our results and give detailed discussions 
on the same. We give our concluding remarks 
in section IV. 



Here /3 = and Q{t) is a time-dependent 
order parameter and is given by 

Q{t) ^ J J dndrip{n,0)p{r^2,t)w{\rl~r^2\) 
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w{r) is the overlap function that is unity inside 
a region of size a and zero otherwise, where a 
is taken on the order of particle diameter. In 
our studies we choose a = 0.40 for all the sys- 
tems with different polydispersity. Q{t) mea- 
sures the number of particles that in a time 
t has either remained within a distance a of 
their original position (when i = j ) or were 
replaced by another particle (when i ^ j ). 
We can separate Q into self and distinct parts, 
Q{t) = Qs{t) + Qoit). The self part corre- 
sponds to terms with i ~ j, Qs{t) = 
ri{0) — rj{t) I). The distinct part is given by 
Qnit) = E r^{0) - Tjit) I)). The sus- 

ceptibility Xi{^) can then be decomposed into 
self, distinct and cross terms Q, 

X4(i)=xf(i)+xf(i)+xr(i) (3) 

where, 

Xlit) (X {Ql{t)) - {Qs{tr (4) 



xf (<) (X {Ql{t)) - {Qn{t))' 



(5) 



II. THEORY AND COMPUTATIONAL 
METHODS 



A. Four-point susceptibility 

The two-point, two-time, fourth-order den- 
sity correlation function 0, 01 is defined as 

g4{ri,f2,t) = {p{fl,0)p{fi,t)p{r2,0)p{ri,t)) 

~{p{rl,0)p{rl,t)){p{r-2,0)p{r-2, t)) 

The volume integral of gi(ri,r2,t) gives the 
four-point susceptibility X4(0i 

= ]^ / drldrip{ri,0)p{fi,t)g4{ri,r2,t) 

(1) 

It has been shown that Xi (t) can be written as 



and, 







X4{t)^^^[{Q'{t))-{Q{tr] (2) 



xfW « {Qs{t)Qi 



it))~{Q sit)) {Qoit)) (6) 

As has been found in previous studies 0], we 
find that for our model system also the ma- 
jor contribution to Xii^) comes from xf (^) and 
hence in this paper we have presented results 
only for xKO- 

The definition of Xii^) in Eq. (1) is in terms 
of spontaneous fluctuations of local dynamics. 
Berthier et al have used fluctuation dissipation 
theorem (FDT) and have defined four-point 
susceptibilities in terms of the response of the 
averaged two-time dynamical correlators to an 
infinitesimal perturbing field [1, [T^ , 



Xxit) 



d{\ fit) 
dx 



(7) 



where x = Tor p. Here d{\ fit) \r) is a standard 
two-time correlator and /(r, t) = 0(r, t)0(r, 0). 
For instance when the observable O is the ex- 
cess density p(r-, t) — ppQ, then d{\ fit) ],-) 
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FIG. 1: The time dependence of the four-point sus- 
ceptibility (X4(i)) a-t four different temperatures 
and at two different polydispersity, 5 = O.fO (thick 
lines) and S — 0.20 (dashed lines). From the bot- 
tom to the top, temperature decreases, xf (*) grows 
for both the systems as T decreases but there is a 
more pronounced growth at lower polydispersity. 



is the intermediate scattering function. Here 
I /(t) |r= y J drf{r, t) is the spatial average 
over a large but finite volume, V . The above 
defintion of Xx{i) provides a very valuable ex- 
perimental tool to measure the dynamic length 
scales in glass systems as shown in [J]. 



B. Bond-orientational order 

The average microscopic structure of liquids 
is usually described by the radial distribution 
function or the structure factor, which essen- 
tially measures only the density-density corre- 
lation function. However, bond-orientational 
order parameters (BOP) introduced by Stcin- 
hardt et al [H, O, |T|] gives a better quantifi- 
cation of the local structure as they capture 
the symmetry of bond orientations. BOP is 
described in terms of combinations of spher- 
ical harmonic functions. Consider a system 
of A'' particles. First, one defines a set of 
"bonds 'which are defined as the vectors con- 
necting neighboring particles. All particles j 
within a cutoff distance tq of particle i arc de- 
fined as neighbors of particle i. Here rg is cho- 
sen to be equal to the distance to the first mini- 
mum of the radial distribution function (RDF) . 
The local order parameters associated with a 
bond r are the set of numbers 

Qirnir)=Yim{0{r),cl)ir)) (8) 



where 9{r) and (f>(r) are the polar and az- 
imuthal angles of the bond with respect to 
an arbitrary but fixed reference frame and 
YimiOiT), 4'{i')) E^rc the spherical harmonic func- 
tions, ft is useful to consider only the even-1 
spherical harmonics, which are invariant under 
inversion. Global bond order parameters can 
be calculated by averaging over all the bonds 
in the system, 

bonds 

Since QimS for a given I depends on the ro- 
tations of the reference frame, it is important 
to consider the rotationally invariant combina- 
tions such as 

Q'^talTT^'^'-''^" ^^^^ 

and 

mi+m2+m3=0 

Wl^ £ ( n^AruAn^s (11) 

mi ,m2 ,m3 

Qi and Wi are the second and third order in- 
variants, respectively. The coefficients (...) are 
Wigner 3j symbols. One also defines a normal- 
ized quantity, 

which for a given I is independent of the mag- 
nitudes of the Qim- The four bond order pa- 
rameters (54, Qq, W4 and Wq are sufficient to 
identify different crystal structures. The typical 
values of these for different crystal structures 
are given in [l^. For a liquid the global values 
of all these four quantities are zero as there is 
no long range order. Note that in clusters with 
cubic symmetry non-zero averages occur only 
for I > 4 whereas non-zero averages occur only 
at Z = 6 and 10 for icosahedral cluster. 



C. System and simulation details 

Micro canonical (NVE) ensemble MD simu- 
lations are carried out in three dimensions on 
a system of = 864 particles of mean diame- 
ter a with polydispersity in both size and mass. 
The polydispersity in size is introduced by ran- 
dom sampling from the Gaussian distribution 
of particle diameters a, 

P{a) = -^exp[-\{'^f] (13) 
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FIG. 2: (a) The value of the peak height of xf (i) is 
plotted as a function of T for S = 0.10 (filled circle), 
S = 0.15 (star) and S = 0.20 (filled triangle) sys- 
tems. The figure shows the suppression of the rate 
of growth of dynamic heterogeneity with S. (b) 
The time at which xf (*) peaks, t* is plotted as a 
function of T for S = 0.10 (SI) and S = 0.20 (S2) 
systems (filled circles and triangles, respectively). 
t* is similar to the a- relaxation time (open circles 
and triangles, respectively). The latter is obtained 
by doing KWW fit to Fs{kmax,t) where kmax cor- 
responds to the first peak in the static structure 
factor. 



FIG. 3: (a) xf (^*) 8- function of the distance 
from the MCT critical temperature Tc for S = 0.10 
and 5 = 0.20 systems. Tc values are 0.42 and 0.39 
for 5* = 0.10 and S — 0.20 systems, respectively. (b) 
a- relaxation time, Tq as a function of the distance 
from Tc 



titles in this study arc given in reduced units 
(length in units of cr, temperature in units of 
1^ and time in units of t = (■^^^)^). 



The standard deviation 6 of the distribution di- 
vided by its mean a gives a dimcnsionless pa- 
rameter, the polydispersity index S — =. The 
mass nii of particle i is scaled by its diame- 
ter rrii = rf%{^)^. We have chosen rn = 1.0. 
The simulations are carried out at different val- 
ues of the polydispersity index, S but at fixed 
volume fraction, cj) = 0.52. The interactions 
between the particles are given by the shifted- 
force Lennard- Jones (LJ) 12-6 potential 

U:, = ^e^Ai^r - i^f] (14) 

where i and j represent any two particles and 
aij = ( °''+°"j The LJ interaction parameter 
e.y is assumed to be the same for all particle 
pairs and set equal to unity. The particles are 
enclosed in a cubic box and periodic boundary 
conditions are used. The cutoff radius Vc is cho- 
sen to be 2.50". A time step of 0.001 is employed 
for T > 1.0 and 0.002 for T < 1.0. AU quan- 



III. RESULTS AND DISCUSSION 

The main objective of our study is to demon- 
strate the effect of polydispersity and hence of 
the local structure on the growth of dynamic 
heterogeneities. By varying polydispersity we 
can Hune'the local structure and hence study 
its effects on the dynamic heterogeneity. As 
polydispersity is increased, the local structure 
is progressively destroyed. Hence the blocking 
of the particles in the cages of the neighboring 
particles (as required for the mode coupling the- 
ory of dynamic transition [l7| ) becomes ineffec- 
tive at higher polydispersity. We find that this 
has a pronounced effect on the development of 
dynamic heterogeneities as well. In this section 
we systemically present our results and show 
that the local structure plays a very important 
role in determining the dynamics in supercooled 
liquids near glass transition. 
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FIG. 4: Angell-like fragility plot at different S. 
The thick lines are VFT fit to the diffusivity data, 
D = Do cap ( y^^^ ). The reference temperature Tr 

is chosen such that D{Tr) = 4.5 x 10"^. The VFT 
extrapolation is used to locate Ty. The plot shows 
that fragility decreases with S. Strength parame- 
ter m (where m = [31) obtained from VFT fit 
have the values 7.78° 8.54 and 15.94 for S = 0.10, 
5 = 0.15 an d 5* = 0.20 systems, respectively. 



A. Suppression of the rate of growth of 
dynamic correlations by polydispersity 

The four-point susceptibility xf (*) obtained 
from Eq. (2) is shown in FIG. [I| for S ^ 0.10 
and S ~ 0.20 for a few temperatures. From 
Eqs. (1) and (2) we see that X4{t) becomes 
larger when the dynamic fluctuations become 
increasingly spatially correlated. Since X4(0 
the volume integral of the four-point correlator 
gi(r,t), it is directly related to the number of 
correlated particles. As temperature is lowered 
xf (i) grows for both the systems but the rate 
of growth decreases with polydispersity. This 
more clearly seen in FIG. [DJa) where the peak 
height of X4(^) (which we label as X4(^*)) is 
plotted against temperature for different values 
of polydispersity. FIG. OJa) shows the suppres- 
sion of the rate of growth of dynamical het- 
erogeneity by polydispersity. In FIG. [3]Ja) we 
plot the peak height X4('^*) ^ function of 
the distance from the MOT critical tempera- 
ture Tc for S = 0.10 and S = 0.20 systems. 
The rescaled plot also shows a suppression in 
the rate of growth of xf (^*) foi' S ~ 0.20 sys- 
tem as compared to S" = 0.10 system. By fit- 



ting to the expression X4(^*) ^ t^" i we 
get the values of the exponent 7^ as 1.507 and 
1.188 for S = 0.10 and S = 0.20 systems, 
respectively. The exponent 7^^. thus seems to 
change with polydispersity. The suppression of 
the rate of growth of dynamical heterogeneity 
by polydispersity leads to the dynamic cross- 
overs observed in the values of the stretch expo- 
nent, (3 and the non-Gaussian parameter, 02 (t) 
between S = 0.10 and S = 0.20 systems as 
shown earlier [ll| and the cross-over behavior 
seen in the exponent z^- that quantifies the de- 
viation from the prediction of Stokes-Einstein 
relation(See Section III. B). Increasing polydis- 
persity at fixed volume fraction decreases the 
fragility of the system (See FIG. Fragility 
measures the rapidity with which the system 
approaches glass transition. Hence decrease in 
the the rate of growth of dynamic heterogeneity 
with polydispersity is consistent with the de- 
crease in the fragility. This is further explained 
in section III.C. 

In FIG. [UJb) we show the time at which 
xfit) peaks, t* versus temperature for S = 
0.10 and S = 0.20. Also shown are the 
a-relaxation times, Tq obtained by doing 
Kohlrausch- William- Watts (KWW) fit to the 
self part of intermediate scattering function, 
Fs{k,t). The plot shows that the dynamics is 
maximally correlated on time scales of the or- 
der of the alpha relaxation time. In FIG. ^h), 
Ta is plotted as a function of the distance from 
the MCT critical temperature, Tc- The rescaled 
plot shows that the rate of growth of Tq, de- 
creases with polydispersity. 



B. Breakdown of Stokes-Einstein Relation 

In this section we discuss the breakdown of 
Stokes-Einstein {SE) relation and its connec- 
tion to the rate of growth of dynamic hetero- 
geneity in the system as temperature is lowered. 
The SE relation is based on treating the liquid 
as a continuum and is given by. 



Here C is a constant that depends on the 
boundary conditions (stick or slip) and rj is 
the viscosity. If the Stokes-Einstein relation 
is strictly valid, then a plot of ln{Di/Dj) ver- 
sus ln(aj/ai) would be a straight line with unit 
slope. Here i and j are indices for solute and 
solvent, respectively. In FIG. [5] we show this 
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FIG. 5: The plot of ln{Di/Dj) versus ln{aj /(Ji) for 
5 = O.lO(circles) and S = 0.20 (triangles). Here 
the subscript i denotes the smaller particle. If the 
SE relation is valid, this plot would be of unit slope 
(dashed line) . The plot shows that at high tempera- 
ture the deviation from the Stokes-Einstein predic- 
tion is higher for S — 0.20 system, but the scenario 
reverses at low temperature where S = 0.10 system 
shows a stronger deviation due to the faster growth 
of dynamic heterogeneity. 



plot for S* = 0.10 and S = 0.20 systems. Both 
systems show deviation from the SE prediction 
even at high temperatures due to the intrinsic 
heterogeneity in the system, with the deviation 
being more pronounced for S = 0.20 system. 
The SE relation has been shown to be not valid 
for the diffusion of small solutes in a solvent 
of bigger particles fT§| . There is an anomalous 
enhancement of the self-diffusion over the SE 
value for small solutes which can be described 
by a power law. 



(16) 



Hence the exponent z^j quantifies the deviation 
from SE relation. It is unity in the SE limit 
and usually larger than unity in supercooled liq- 
uid. FIG. [6] shows that deviates significantly 
from unity for polydisperse liquids, particularly 
at low temperatures. Interestingly, it is larger 
than unity even at high temperature because 
of the heterogeneity in the size and mass. This 
deviation of the slope from unity at high tem- 
perature can be a combination of two different 
effects. The first one is the mass which is not 
present in SE relation but has been reported 
earlier in simulations [l9l and mode coupling 
theory (MOT) studies [201. The studies predict 
a weak power law mass dependence of diffusion. 
The second effect is that of size which has also 
been obtained in experiments and simulations 
in and MCT studies [li . When the size of one 



of the particles is 1.5 — 15 times smaller than the 
other it shows an anomalous rise in diffusion. 
This enhanced diffusion has been explained in 
terms of microviscosity effect. The MCT stud- 
ies explain the microscopic origin of the size 
effect in terms of the difference in relaxation 
timescales of the two particles which leads to a 
decoupling in the dynamics [3| ■ (Note that for 
size dependent studies the small particle was a 
tracer and for mass dependent studies the heav- 
ier particle was a tracer. In the present study 
the systems arc intrinsically heterogeneous.) 

As temperature is lowered the deviation from 
the prediction of SE relation becomes more 
pronounced and one observes a crossover in the 
value of Za between the 5* = 0.10 and S — 0.20 
systems; the values are much higher for S = 
0.10 system than S = 0.20 at low temperatures 
[see FIG. [B]. This again shows that the rate 
of growth of dynamic heterogeneity is faster in 
S — 0.10 system than in S* = 0.20 system. 
The faster rate of growth of dynamic correla- 
tions leads to similar temperature-dependent 
crossovers between S — 0.10 and S ~ 0.20 sys- 
tems in the values of the stretch exponent (/3) 
and the non-Gaussian parameter [a2{t))^ both 
of which contain implicit information on dy- 
namic heterogeneity [ll|. These studies show 
that there is a strong correlation between the 
growing dynamic heterogeneity in the system 
and the breakdown of SE relation as the for- 
mer renders the continuum description of liquid 
invalid as required for the SE relation. 



Fragility and the growth of dynamic 
correlations 



It has been shown that the dynamics of 
fragile liquids are more spatially heterogeneous 
than that of strong liquids [2^. Increas- 
ing polydispersity at fixed volume fraction de- 
creases the fragility of the system as shown in 
FIG. m The decrease of the rate of growth 
of X4(^) with polydispersity supports the pre- 
viously observed correlation between fragility 
and dynamic heterogeneity. The intrinsic het- 
erogeneity of the system (as measured by the 
distribution of particle masses and sizes) in- 
creases with polydispersity. Hence we have the 
interesting scenario in which increasing polydis- 
persity leading to a more homogeneous dynam- 
ics even though the system becomes completely 
amorphous at higher values of polydispersity. It 
has been shown that polydispersity has a pro- 
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FIG. 6: The values of thepower law exponent, Za 
obtained by fitting Eq. f6, are plotted against T, 
for S = 0.10 and S = 0.20 systems. If the SE 
relation is strictly valid then Za = 1- Deviation 
from unity shows the breakdown of SE relation. 
The figure clearly shows that at high temperatures 
intrinsic heterogeneity causes a larger breakdown 
in S = 0.20 system, whereas at low temperatures 
the faster growth of dynamic heterogeneity leads to 
a stronger breakdown in S = 0.10 system. 



nounced effect on potential energy surface [23[ . 
As polydispcrsity is increased from zero, the 
characteristics of the potential energy minima 
change from that of crystalline to that of amor- 
phous. The latter is known to have low curva- 
ture and small barriers along some coordinates 
[13i [13 ■ This observation is also consistent 
from the perspective of the inherent structure 
formalism, according to which the potential en- 
ergy landscape of a fragile liquid is very hetero- 
geneous which in turn leads to heterogeneous 
dynamics whereas the landscape of strong liq- 
uids consist of a single mega basin [2^ . Hence 
from a potential energy landscape perspective, 
increasing polydispcrsity leads to a smoothen- 
ing of the landscape that in turn leads to the 
facilitation of dynamics as well as decrease of 
fragility. In Section III. D, we try to under- 
stand how polydispcrsity suppresses the rate of 
growth of dynamic correlations and in particu- 
lar, whether the loss of structure upon increas- 
ing polydispcrsity has any role to play in this. 
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FIG. 7: The calculated static structure factor, S{k) 
is plotted against wave number k for S = 0.10 
(thick lines) and S = 0.20 (dashed lines) systems 
at a few temperatures. Structural correlation is 
weaker in S = 0.20 system than in S = 0.10 system 
and shows no appreciable change with temperature. 



D. Local structure and the growth of 
dynamic heterogeneities 



We plot tlie static structure factor, S{k) for 
S = 0.10 and S = 0.20 systems in FIG. [3 The 
plot shows that increasing polydispcrsity de- 
stroys the local structure in the system as the 
system becomes more amorphous. The peak 
height of S{k) is highly suppressed in S = 0.20 
system as compared to S = 0.10 system and 
does not show any appreciable growth upon 
lowering of T. FIG. [5] shows the peak height 
value of RDF, g(rmax) as a function of temper- 
ature. At S" = 0.10, the peak height shows con- 
siderable enhancement upon lowering of tem- 
perature whereas at S = 0.20 there is no re- 
markable change in the value of g{rmax) with 
temperature. 

As mentioned in Section II. B, the bond ori- 
entational order parameters give a better quan- 
tification of the local structural arrangement. 
Frank [23| proposed that atoms might form 
icosahedral clusters in liquids since the lowest 
energy state of a 13-atom cluster interacting via 
Lennard-Jones potential is an icosahedron (and 
not fee). But icosahedra cannot tile space in 
3-dimensions due to its 5-fold symmetry and 
hence do not satisfy the global structural stabil- 
ity criterion. This geometrical frustration could 
be an important factor that contributes to the 
stability of glassy state [11]. Steinhardt et al 
have shown that there is a long range ori- 
entational icosahedral order in supercooled liq- 
uids. It has been shown that the large size dis- 
parities at higher values of polydispcrsity would 
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FIG. 8: Peak height value of the radial distribution 
function, g{rmax), when plotted against separation 
r, is plotted as a function of T for S = 0.10 (cir- 
cles), 5 = 0.15 (stars) and S = 0.20 (triangles). 
S = 0.20 system does not show any remarkable 
change in the value of gifmax) upon lowering of 
T. On the other hand, S = 0.10 system shows a 
sudden increase of spatial correlations for T < 0.8. 
Comparison between FIG. [8] and FIG. EJa) shows 
that local structure plays a crucial role in the build 
up of dynamic correlations. 



inhibit any icosahedral cluster formation. How- 
ever, since at low/moderate polydispersity the 
peak height of g{r„iax) shows a pronounced 
growth as temperature is lowered, one can ask 
whether this is due to the formation of icosahe- 
dral clusters that grows with decrease in tem- 
perature. 

We look for the local values of BOP in order 
to understand whether local orientational order 
plays any role in the growth of dynamic het- 
erogeneities. The icosahedral order, if present, 
would be picked by the BOP corresponding to 
I = 6, Qe- The local values of Qe are plotted 
in FIG. O To get the local values the spherical 
harmonics corresponding to / = 6 are summed 
over the nearest neighbor bonds only. The fig- 
ure shows that there is a pronounced icosahe- 
dral orientational order at the local level. This 
local icosahedral order shows considerable en- 
hancement at lower polydispersity as tempera- 
ture is lowered (for T < 0.80) and also decreases 
with polydispersity. In the inset of FIG. [9] we 
plot Qq as a function of the distance from the 
MCT critical temperature Tc which shows that 
as temperature is lowered, Qe increases much 
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FIG. 9: The calculated local values of bond or- 
der parameter, Qa are plotted as a function of 
T for S = 0.10 (circles), S = 0.15 (stars) and 
S = 0.20 (triangles). The plot shows that at the 
local level there is significant orientational order 
that increases with decrease in temperature and de- 
creases with polydispersity. [Inset: Qe as a function 
of the distance from Tc-] 



more sharply at lower values of S. 

In FIG. [To] we plot the global values of Qq for 
different 5 as a function of temperature. The 
averages over bonds are evaluated by summing 
over all bonds lying within a sphere of radius 
2.4 units. Nine such spheres are considered 
whose centers lie at different locations of the 
simulation box. We repeat this averaging for 
several different snapshots obtained from sim- 
ulation. It is evident from FIG. [10] that poly- 
dispersity suppresses long range orientational 
order. Even at moderate polydispersity, there 
is no pronounced growth of long range icosahe- 
dral order upon supercooling. 

Our results indicate that increasing polydis- 
persity destroys both the local structure and 
the local orientational order. The four-point 
susceptibility xf (i) measures the susceptibility 
arising from the number of localized particles 
and is a measure of the dynamic heterogeneity 
in the system. Thus the dynamic heterogeneity 
is associated with the temporary localization 
of particles by their neighbors. Since the lo- 
cal structure is destroyed, at higher values of 
polydispersity it is not possible to have such a 
caging effect. As a consequence particle motion 
gets decorrelated over much shorter time scales. 
This is best seen by plotting the van Hove cor- 
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FIG. 10: The calculated values of the global bond 
order parameter, Qg are plotted as a function of T 
for S = 0.10 (circles), S = 0.15 (stars) and S = 0.20 
(triangles). The plot shows that there is no appre- 
ciable long range orientational order developing in 
the supercooled state. 




FIG. 11; The calculated van Hove self-correlation 
function, Ga{r,t), plotted against position r at dif- 
ferent times (indicated on the figure) for 5 = 0.10 
(thick lines) and S — 0.20 (dashed lines) systems at 
T = 0.45. The figure shows the faster decay of den- 
sity correlations for 5* = 0.20 system as compared 
to 5 = 0.10 system. 



E. Polydispersity and Glass Forming 
Ability (GFA) 



Wc find that the present system of LJ parti- 
cles of varying size and mass crystallizes when 
polydispersity is less than 7.5%. Thus, our sys- 
tem with 10% polydispersity can be regarded as 
the system on the lower side of polydispersity 
that could be made to avoid crystallization and 
remain liquid within our MD simulation time 
range. Interestingly, we find that this is also 
the system that shows glassy behavior at the 
highest temperature. When we increase poly- 
dispersity beyond 10%, we need to lower the 
temperature to capture the onset of slow glassy 
dynamics. This can be seen from FIG. (Hb) 
which shows that the rate of growth of Tq de- 
creases with S. The new aspect revealed in 
the present work is the correlation between the 
GFA and the rate of growth of the dynamic 
heterogeneity — sharper the growth, larger the 
GFA. 

Given that a polydisperse liquid with low 
polydispersity {S < 0.05) crystallizes easily, the 
loss of local structure at large S (> 0.20) and 
the concomitant difficulty of glass formation at 
large S imply a rather narrow range of S for 
polydisperse systems to act as good glass for- 
mers. This means that only at moderate poly- 
dispersity the system has a high GFA. The 
GFA decreases with polydispersity beyond a 
value of S. Further insight can be gained from 
the study of inherent structures. We find that 
the ruggedness of potential energy landscape 
decreases with 5', which is consistent with de- 
crease of the GFA as well as fragility with S. 
It is important to note that network glass for- 
mers like Silica which is a strong liquid in An- 
gells fragile/strong classification, exhibits high 
glass forming ability due to trapping by de- 
fects. This apparently contradicts the decrease 
of GFA with polydispersity. The latter appears 
to be a hallmark of polydisperse systems. We 
shall address these issues in detail elsewhere. 



IV. 



CONCLUDING REMARKS 



relation function [see FIG.s [TT] and [12] . Thus 
the loss of local structure due to polydispersity 
suppresses the growth of dynamic heterogeneity 
in the system. 



The hypothesis that structure determines dy- 
namics has been termed by Harowell as the cen- 
tral dogma of glass science^^. This dogma is 
validated in the Mode Coupling Theory. The 
Adams- Gibbs theory, however, gives larger em- 
phasis on the emergence of a dynamical cor- 
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FIG. 12: The calculated van Hove distinct- 
correlation function, Gd{r,t) is plotted against po- 
sition r for 5 = 0.10 (upper panel) and 5 — 0.20 
(lower panel) systems, at T = 0.45 depicting the 
faster decay of inter-particle correlations at higher 
polydispersity. 



relation length as the source of slow dynamics 
which does not seem to depend too sensitively 
on the structure formation. This can be un- 
derstood from the relative inscnsitivity of the 



structure to temperature. In the present work, 
we have varied polydispersity that allows large 
variation of the local structure and find that the 
local structure indeed plays an important role 
in the development of dynamic correlations and 
the slow dynamics near glass transition in a su- 
percooled polydisperse liquid. Increasing poly- 
dispersity at constant volume fraction leads to 
a suppression of the rate of growth of dynamic 
heterogeneity in the system, which can be at- 
tributed to the loss of local structure with poly- 
dispersity. At moderate polydispersity, there is 
a faster growth of structural correlations as the 
temperature is lowered, which leads to a corre- 
sponding faster growth of dynamic heterogene- 
ity. At higher polydispersity, structural correla- 
tions arc weak and do not show any significant 
change with temperature and correspondingly, 
the rate of growth of dynamic correlations is 
also less. We also find that there is a pro- 
nounced local icosahedral order which increases 
with cooling and decreases with polydispersity. 
No significant long range icosahedral order is 
found either in the equilibrium or supercooled 
liquid. 

An important outcome of the present work is 
the hitherto unknown correlation between poly- 
dispersity and glass forming ability. This cor- 
relation deserves further study. 
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